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Abstract 

We study both massless scalar and Yang-Mills field theories in the deep infrared in presence of 
a simple boundary. We can show, with the help of the recent scenario emerging from studies on 
the propagators, that the presence of a mass gap makes the Casimir contribution exponentially 
small as should be expected. The existence of a trivial infrared fixed point shared by both theories 
makes the computation as simple as the free particle case. 
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I. INTRODUCTION 



Vacuum in quantum field theory is known to not be inert. The most striking evidence 
of this was put forward on 1945 by Casimir [l| providing an example of a macroscopic 
effect produced by quantum fluctuations. Due to the smallness of the involved elements, 
an experimental proof of the Casimir effect has been obtained quite recently (see [2j for a 
review) . 

Casimir effect is an example of a non-trivial behavior of a quantum field theory due to 
the presence of boundary conditions that somehow modify the behavior of the free space 
case. This kind of question is interesting per se. But while in the ultraviolet regime we 
have techniques able to cope with this kind of problems, in the opposite limit, the infrared 
regime, such techniques were generally lacking. It is important to note that the behavior at 
different energy regimes of a quantum field theory could result quite different depending on 
the structure of vacuum. So, we know that at high energies, QCD behaves as an almost free 
theory and asymptotic states of the Yang-Mills field are massless gluons. On the other side, 
at lower energies, a Yang-Mills theory displays a mass gap yielding a very different behavior 
in this case. E.g. this could be inherent to a non-trivial non-perturbative vacuum seen as 
an instanton liquid 



lad an important rebirth in 

n 

5[ and by the improvement 



Studies of quantum field theories in the infrared limit have 
this last two decades by the use of Dyson-Schwinger equations [4J, 
in computing resources that made possible to analyze Yang-Mills theory on larger lattices. 
The theoretical proposal to solve the infinite hierarchy of Dyson-Schwinger equations was a 
truncation scheme that, for the Landau gauge, produced a gluon propagator going to zero 
and a ghost propagator running to infinity faster than the free case with momenta going 
to zero. Such a scenario was previously devised by Daniel Zwanziger jf| providing also a 
criterion for color confinement then dubbed Gribov-Zwanziger criterion. Initial lattice com- 
putations seemed to support these conclusions even if no bending toward zero of the gluon 
propagator was ever observed at lower momenta. People generally thought that was just a 
matter of volumes and, increasing computational resources, things should have changed. 

The breakthrough come out on 2007. At the Lattice 2007 Conference in Regensburg, 
three groups presented their results with huge volumes arriving to such a significant value 
as (27/m) 4 [7H9|. The shocking result was that the scenario devised since then, generally 
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accepted as correct, was not describing the situation seen on the lattice: The gluon propa- 
gator was reaching a plateau at lower momenta with a finite non-zero value in zero and the 
ghost propagator was behaving as that of a free massless particle. All in all, the running 
coupling was seen to bend clearly toward zero without evidence of a non-trivial fixed point 
as was expected instead. 

In the eighties, a classical paper by Cornwall 10] showed that indeed the gluon propagator 
showed a dynamical mass, dependent on momenta, that in the zero limit reaches a finite 
constant value. With the emerging of techniques to solve Dyson-Schwinger equations in the 
nineties, it was straightforward to try to solve them numerically. This numerical approaches 
showed that Cornwall view is indeed correct 11( but this paper displayed also all the scenario 



seen since 2007 on lattice computations. So, a research line developed producing an in-depth 
theoretical and numerical analysis of Dyson-Schwinger equations supporting the current view 
12-2J. 

The idea in this paper is to start from this situation, also supported by theoretical 
arguments, to analyze the behavior of the vacuum in the infrared limit. The best way 
to see a non-trivial behavior is through the analysis of the deep infrared limit with simple 
boundary conditions. As we will see, this can be accomplished yielding the result, somewhat 
expected with a mass gap, that the Casimir contribution is exponentially damped both for 
scalar and Yang-Mills theories. 

The paper is so structured. In sec Jill we discuss the theoretical approach to treat quantum 
field theories in the infrared limit. In sec JIIII we compare our results with lattice compu- 
tations and numerical Dyson-Schwinger equations showing how these strongly support our 
conclusions and providing a theoretical framework to perform our computations. In sec JIVI 
we present our main results on the Casimir contribution. Finally, in secJV] we give our 
conclusions. 



II. INFRARED LIMIT 

Infrared physics can be studied through perturbation techniques as already shown in 



22 



-|24j much in the same way as it is seen in a weak coupling case. The results appear 
to be quite enlightening producing explicit analytic solutions both for the quantum and 
the classical cases in the limit of a bare coupling taken to go formally to infinity. This 
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approach can be seen as a substantially improved version of the approach devised at the 



end of seventies by Carl Bender and others [25 



26l |. The relevant result useful for our aims 



is that both the scalar and Yang-Mills theories display a trivial infrared fixed point making 
really easy to determine the effective potential in this limit. 



A. Scalar field theory 



For our aims it is enough to consider the following generating functional of a scalar field 



Z[j] = N J [d(f)]expli J d 4 x 



(1) 



One can formally rescale space-time coordinates as x — > \/Xx and the functional can be 
rewritten as 



Z[j] =N J[d</>] expli J d 4 x 



—my - V+ - 



(2) 



We note here that j is an arbitrary function and can be rescaled as j/X — > j. So, at this 
stage, we can perform an expansion 



n=0 



(3) 



producing the following terms into the action 

S = J d 4 x 
S\ = J d A x 
So = f d 4 x 



(4) 
(5) 
(6) 



From these, we are a step away from a proof of triviality of this theory in the infrared 
limit. Indeed, we just note that this functional demands, undoing the rescaling, to solve the 
equation 

«9 2 0o + = j (7) 
in the limit A — > oo. Using this equation, we get a functional 

Z\j] » e < / rf M3W0M*S+i*>] /[^e'j/^IIW-Wl (8) 

where we note that the relation O — 4>o\j\ is hidden inside the functional integral. To work 
out this relation, we note that we can solve eq.flT]) by taking the following approximation 



being 



0o = yU / d x'A(x — x')j(x') + 5(f) 



d 2 A(x - x') + \[A{x - x')f = -5\x - x'). 



(9) 



(10) 



with fi an arbitrary constant with the dimension of energy. This constant will define the 
mass spectrum of the theory as we will see below. One gets 



ti\ J d 4 x'A(x - x') j/i J d 4 x"[A(x' - x")fj(x") 



n3 



d A x"A{x' - x")j{x") 



+ ... 



(11) 



and higher order corrections can be obtained by iteration. This appears a functional ex- 
pansion in the current j that is consistent with the large coupling approximation. The 
propagator can be computed exactly [22|, |29| and can be expressed in the following form 



Hp) = E 



being 



with K(i) = Jfj 



B 



(2n 



n=0 p 2 -m 2 n + it 



71-2 (_l)n e -(n+i)7r 



a/ 1+sin 2 e 



' K 2 {i) i + e -( 2n + 1 ^ 
1.3111028777 and 



(12) 



(13) 



rrir. 



1 

n H — 



7T 



(14) 



1) K{%) \2 j 

the mass spectrum of the theory, the one of a harmonic oscillator. This can be easily seen 
as the Fourier transform in time at p = of this propagator has the form (0(0, t)<fr(0, t')) = 
J2n°=o C n e~ m "( t ~ t '' ) showing that, in the infrared limit, the theory develops a mass gap. This 
is the key result of our analysis that also completes our proof of triviality of this scalar field 
theory in four dimensions in the infrared limit. The final form of the generating functional 
at the leading order can be written down as 



Z [j] = iVexp 



- d xd yj(x)A(x - y)j(y) 



(15) 



having the expected Gaussian form. We observe that the free excitation we have found in 
this limit entails a subset of excited states with the spectrum given by eq. (I14j) . For the 
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sake of completeness, we give here the next to leading order correction to this generating 
functional 



Z\j\ « Z [j] J [#] exp {I J d 4 x [i(50) 2 - (/ rf 4 XiA(x - x 1 )j(x 1 )^j 2 



(16) 



B. Yang-Mills theory 



Firstly, we review an approach devised in the eighties 30|, |31[ that gives a clear under- 
standing of the stumbling block arisen in the studies of infrared Yang-Mills theory. Quantum 
Yang-Mills theory can be stated in its simplest form through the following generating integral 



i J d 4 x (j;^ a + r c a + cv) 



^TrF 2 + — (<9 ■ A) 2 + {c a d^c a + g^f ab %A b »c c 



(17) 



with the field strength given by F£ u = d IJ ,A^ — d l ,A^+gf abc A b A^ and A^ the vector potential. 
Now, we can proceed as for the scalar field. In this case we just rescale x — > yNgx being N 
the number of colors. Then, in order to understand the behavior of Yang-Mills theory in the 
infrared limit, we need a way to manage the classical equations of motion (with standard 
notation) 

= -3 a v (18) 
being = <9 M — igt a A a ^. Indeed, we can write down the solution in the form 

A;(x) =Aj d'yD ab v (x - y)f b (y) + O (l/VNg) (19) 

being A a constant having the dimension of a mass and this asymptotic approximation will 
be clearer in a moment. In this way we recover a Gaussian approximation at the leading 
order and a trivial fixed point at infrared. But at this stage we have only a conjecture unless 
we are able to put forward the exact form of the gluon propagator in the infrared limit, 
the stumbling block we talked about above. This question can be approached through a 
mapping theorem proved recently 



32 



33j that can be stated in the following form: 



Theorem 1 (Mapping) An extremum of the action 

S = Id'x I(50) 2 -^0 4 



(20) 
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is also an extremum of the SU(N) Yang-Mills Lagrangian when one properly chooses with 
some components being zero and all others being equal, and A = Ng 2 , being g the coupling 
constant of the Yang-Mills field, when only time dependence is retained. In the most general 
case the following mapping holds 



A;(x)=r ] ;<j ) (x)+0(l/VNg) (21) 
being r]® constant, that becomes exact for the Lorenz gauge. 



The proof of this theorem was completed in Ref. 33] to answer a criticism by Terence Tao. 



Tao finally agreed with this conclusion 34| . Stated otherwise, this theorem determines an 



asymptotic mapping between the solutions of the two classical theories when the couplings 



are taken large enough. An incomplete form of this mapping was already stated in Ref. [35 1 . 
The existence of this mapping grants that we can write down the propagator of Yang-Mills 
theory, e.g. in the Landau gauge, as 

D±{V) = <U - ^] A(p) + O (-L-) (22) 



™ V p 2 J - WN 9/ 

with A(p) given by eq. (112ft . provided we change A — >■ Ng 2 and, for our aims, we dub the 
constant /x as A so that the spectrum of a SU(N) can be written down as 

m n =(n + I) -J- (^-] 4 A. (23) 



2J K(i) V 2 

This result is very easy to prove from the definition of the two-point function 

D$,{x-y) = {TAl{x)Ai(y)). (24) 

Applying the mapping theorem above one gets immediately the behavior of the theory in 
the infrared 



= TifalA(x-y) + 0(l/VNg) (25) 

that produces the result. We note the identity for the Landau gauge as a possible Smilga's 
choice 

V a d = 8* L, - P -f\ (26) 



being 77^ the Minkowski metric. 

For the ghost, we note that the mapping theorem grants its decoupling at the leading 
order. This can be seen immediately from the action of the Yang-Mills theory by direct 
substitution. So, at the infrared fixed point, the propagator of the ghost field is that of a 
free particle and we can write 

G(p) = ^— + 0(l/VNg). (27) 
Finally, for the generating functional one has 

Z\j\ = N'exp^Jd 4 xd'yr a (x)D; b u (x-y)f b (y) 

+ O (-±J\ (28) 



'Ng, 

that takes also a Gaussian form showing that Yang-Mills theory in the infrared limit displays 
a trivial fixed point. This is the limit we are interested in here. We note as the spectrum 
of the theory at the trivial infrared fixed point, given by eq. (|23l) . is made of free massive 
excitations, with a superimposed spectrum of a harmonic oscillator, and so entails a mass 
gap. These can be considered as asymptotic states to start from to build up a perturbation 
theory in the infrared limit. This same conclusion does not hold in QCD as this theory has 
a non-trivial infrared fixed point due to the presence of quarks. 

III. COMPARISON WITH NUMERICAL COMPUTATIONS 

In this section we show how sound is the choice of the propagator describing low-energy 
physics starting from measurements obtained from lattice computations and numerical so- 
lution of Dyson-Schwinger equations. This kind of computations, relatively to the Landau 
gauge, span a two decade long period that has seen its main breakthrough on 2007 with the 
clear evidence that the gluon propagator in the Landau gauge is sitting on a plateau at very 
low-energy, reaching a finite non-zero value at zero momenta [t- 9] as also was proven in 11 1 



for numerical Dyson-Schwinger equations. This means that, in order to show the soundness 
of our results given in the preceding sections, we have to compare with these computations. 

We consider two kind of lattice computations: A set of volumes till 80 4 directly obtained 
with measurements on the lattice and measurements at 128 4 recovered from figure 2 in Q|. 
We are able to show in this way that, increasing the volume, our propagator describes even 

8 



more accurately the one measured on the lattice in the deep infrared. We would like to 
point out that the mass gap is different for these two cases as it depends on the value of 
(3. Then, using numerical Dyson- Schwinger results 13) where no volume problem arises, we 
see that our propagator perfectly matches the numerical solution in the deep infrared and 
deviates from it in the intermediate regime where our approximation is expected to worsen. 
Note that we consider a weak dependence on the gauge group as showed in 36] that is fully 
consistent with our discussion above. 




(3=6 (48)" 
p=6 (64)" 
p=6 (80)" 
P=6.2 (48)" 
p=6.2 (64) 4 
p=6.2 (80) 4 



FIG. 1: Gluon propagator in the Landau gauge for SU(3), 80 4 with a mass gap of mo = 287 MeV 
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FIG. 2: Gluon propagator in the Landau gauge for SU(2), 128 4 with a mass gap of mo = 435 MeV 
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FIG. 3: Gluon propagator in the Landau gauge for SU(3) obtained by numerically solving Dyson- 
Schwinger equations and a mass gap mo = 314 MeV 

This agreement between lattice computations at increasing volume and the perfect match 
for the numerical Dyson- Schwinger equations with our propagator give a strong support to 
the idea that a pure Yang- Mills theory reaches an infrared trivial fixed point. This same 
conclusion cannot be drawn for QCD due to the presence of quarks. 

IV. EFFECTIVE POTENTIAL 

In order to compute the effective potential for the Casimir effect, we consider the simplest 
geometrical setting. We assume infinite plane boundaries for the axes x,y and periodic 
boundary condition on the z coordinate. In our case, it is not difficult to write down the 
expected effective potential having an action made up of a sum of free contributions. 

A. Scalar field theory 

The action has the very simple form in Euclidean metric 
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showing a sum of weighted free propagators. Now, this gives us the effective action 

J oo 

T = --Y,B n logdet(p 2 + m 2 n ) (30) 

A n=0 

and using the standard relation In det = Tr log, with the given boundary conditions one has 
T f dpod Pl dp 2 ~ / 2 2 2 /2n,7r\ 2 \ 

being L the distance between the two slabs, S the surface of the boundary. In order to get 
the Casimir contribution we have to evaluate this expression. Firstly, we note that the sum 
on n z is well-known in thermal field theory as a Matsubara sum j37|. So, one has 

l - £ log (u 2 n + (j^j ^ = ^ + In (l - e~ L ^) + constant (32) 



being u n = ypg +Pi +P2 + m n an d the constant is divergent but independent on L or uj. 
Then, the Casimir term is straightforwardly obtained as [2j 

oo . 

ML) = E B n / (2^3 In (l - c-^) S. (33) 

This is the standard result if there is no mass spectrum and no mass gap, that is m n — > 0. 
In the infrared limit, being m n finite, we have to pursue a different approach. We should 
consider small momenta and a cut-off /i, the same as in the mass spectrum. Then, we are 
able to evaluate the integral. One has 

oo 3 3 

E C (L) w-£ B n c~ Lm ^S « -B e- Lm °f^S (34) 

where the approximation Lm 3> 1 for the mass gap has been used. We can conclude from 
this equation that, in the infrared limit, the Casimir contribution is exponentially small. 
This result should be expected from the simple fact that the appearance of a mass gap 
makes forces short-ranged. 



B. Yang-Mills theory 

Due to the existence of a mapping theorem between scalar and Yang-Mills fields we can 
draw a similar conclusion for a gluonic field. This conclusion can change when either quarks 
are present or a mass gap goes to zero increasing momenta scale recovering asymptotic 
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freedom. This makes interesting a study at finite temperature with given boundaries. So, 
taking at the infrared trivial fixed point the action 

So = ~\ J ^r( P )D<?Mf b (-p) (35) 

being 

D%(x-y) = 6 ab ^ u -^A(p) (36) 

the gluon propagator in the Landau gauge and A(p) the propagator of the scalar field. But 
if we use current conservation the action is immediately reduced to 

So = ~\ J ^f a (p)mjt(-p) (37) 

that maps perfectly on the action of the scalar field provided a multiplying factor d(N 2 — 1) 
that is exactly the number of current products. This means that the result is identical to 
the one of the scalar field apart an inessential numerical factor and we can conclude that 
the Casimir effect is exponential small also in Yang-Mills theory as expected. 



V. CONCLUSIONS 



Using recent analyses, both theoretical and numerical ones, on the behavior of the propa- 
gators of quantum field theories in the infrared limit made quite simple to derive the behavior 
of the corresponding vacuum in presence of a simple boundary. In the deep infrared we must 
expect a quite different behavior for a Yang-Mills theory due to a non-trivial vacuum, gener- 
ally non-perturbatively well-described by a liquid of instantons, due to the appearance of a 
mass gap in the theory. A massive scalar field in the free case shows a Casimir contribution 
exponentially small. Here, we can see the same phenomenon to appear as both a massless 
scalar field theory and Yang-Mills theory acquire mass dynamically displaying a mass gap. 
But the situation and even more better due to the fact that these theories hit a trivial 
infrared fixed point making for us very easy to adopt the approaches for free theories. In 
the near future will be very interesting to extend this analysis to more complex boundaries 
and include quark matter. 
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